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Abstract 

We calculate a single transverse-spin asymmetry for the Drell-Yan lepton-pair's 
angular distribution in perturbative QCD. At leading order in the strong cou- 
pling constant, the asymmetry is expressed in terms of a twist-3 quark- gluon 
correlation function Tp\xi, X2). In our calculation, the same result was ob- 
tained in both light-cone and covariant gauge in QCD, while keeping explicit 
electromagnetic current conservation for the virtual photon that decays into the 
lepton pair. We also present a numerical estimate of the asymmetry and compare 
the result to an existing other prediction. 

PACS: 13.85.Qk,13.88.+e 
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I. INTRODUCTION 



It has been known for some time that cross sections involving a single transverse hadronic 
spin show significant asymmetries relative to the spin direction MM. Because of parity and time- 
reversal invariance, single longitudinal-spin asymmetries vanish for most inclusive observables 
in high energy collisions; and single transverse-spin asymmetries open up a new domain of QCD 
dynamics: physics at twist-3, which is sensitive to the correlations between quark and gluon 
fields. With the active spin program at Brookhaven National Laboratory's relativistic heavy 
ion collider (RHIC), many single transverse-spin asymmetries have been proposed and could 
be observed in the near future || |f| . 

In this paper, we re-examine the single transverse-spin asymmetry in the Drell-Yan lep- 
ton angular distribution. Although the asymmetry was studied in Refs. two different 

analytical expressions were derived. In Ref. || a formal argument, based on gauge vector in- 
dependence, was given in favor of the result of Ref. || , however, uncertainty about the result 
remained fill. This motivates us to calculate the asymmetry in a more transparent way, such 



that color and electromagnetic gauge invariance are manifest. 

For the polarized Drell-Yan process, A(p, st) + B{p') — > J*(Q)[— > ££] + X, we define the 
lepton-pair's angular distribution in the virtual photon's rest frame (Q = 0). As shown in 
Fig. [TJ, we choose the z-axis along the direction of polarized incoming hadron, x-axis along the 
direction of polarization vector st, and the observed lepton momenta as 

= Q (1, sin 9 sin 0, sin 9 cos 0, cos 9) ; 

= Q (I; — sin 6* sin 0, — sin 9 cos 0, — cos#) . (1) 

The Q is the invariant mass of the virtual photon, and its four-momentum = + . 

The Drell-Yan single transverse-spin asymmetry is defined to be the difference of two spin- 
dependent cross sections, (1<j(st) / dQ 2 d£L with opposite directions of polarization, divided by 
their sum, 



.4 



N 



1 / cI<j(st) da(—sx 

2 \dQ 2 dQ 

dAa(sr 
dQ 2 dtt 




(2) 

where the phase space dQ = dcos9d<p. 

In Ref. 0, Hammon, Teryaev, and Schafer (HTS) calculated the Drell-Yan single transverse- 
spin asymmetry at tree-level in perturbative QCD with the light-cone gauge, and obtained 
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where g = y/Ana s is the coupling constant for strong interaction, and the summation over 
quark flavor was suppressed. The q(x) is the normal unpolarized quark distribution. The 

was called a twist-3 soft gluon pole function, with the Tp{x%, X2) 

defined as [HI 
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dyi d V2 

47T 

x (p,s T |^(0)7 + e? Tann F a +fa) V(2/r)b,s r >. (4) 



T£ V) (xi,X 2 ) = J c ixiP + vT+i(x2-xi)p+y; 



The ± components of any vector k are defined in terms of two lightlike vectors n and n as 
k + = k ■ n and k~ = k ■ n. These lightlike vectors are chosen such that n-n—1. In our choice 
of frame, the momentum of the polarized hadron p is up to a mass term proportional to n. 

In Ref. ||, Boer, Mulders, and Teryaev (BMT) also investigated the Drell-Yan single 
transverse-spin asymmetry. A similar but different result was derived. This result was re- 
expressed in the notation of Ref. |7J and was found to be [Q 



a (BMT) _ 
N — 9 



sin 29 sin 
1 + cos 2 e 



I T(x, x) . . 



which differs from A^ TS ^ in Eq. (||) by the term proportional to x dT(x,x)/dx. 

In this paper, we calculate the leading order contributions to the Drell-Yan single transverse- 
spin asymmetry in terms of a generalized QCD factorization theorem ||[lT|. in order to test 
the gauge invariance of our result, we derive the asymmetry in light-cone gauge as well as in 
covariant gauge of QCD while keeping the electromagnetic current conservation of the virtual 
photon explicit. The latter also makes the calculation much more transparent than the one of 
Ref. ||, which uses so-called effective T-odd distribution functions. In our present calculation 
we find that the results are independent of the gauge choices and we confirm the result of Ref. 
II- 

The derivative term in HTS's result in Eq. has been found in single transverse-spin 
asymmetries in other processes, like prompt photon production || and pion production in 
hadron-hadron scattering ||, where the derivative term arises from collinear expansions of the 
partonic parts while calculations were performed in covariant gauge. However, we argue below 



that in terms of the generalized factorization theorem |TT|, there is no derivative term for the 
Drell-Yan asymmetry at the tree-level while it can appear in higher order corrections. 

The rest of this paper is organized as follows. In the next section, we derive the single 
transverse-spin asymmetry in the Drell-Yan lepton angular distribution in perturbative QCD 
with a choice of light-cone gauge. In our calculation, we monitor the electromagnetic current 
conservation of the virtual photon by calculating the strong interaction part and the leptonic 
part separately. In addition, we show explicitly why there is no derivative term. As a cross 
check, we use the method developed in Ref. || to calculate the same Drell-Yan single transverse- 
spin asymmetry in a covariant gauge in Sec. III. An identical result is obtained. Finally, we 
present a numerical estimate of the asymmetry and our conclusions in Sec. IV. 



II. ASYMMETRY CALCULATED IN QCD LIGHT-CONE GAUGE 

In order to calculate the Drell-Yan single transverse-spin asymmetry, A # in Eq. (|2]) , we need 
to calculate both the unpolarized cross section, da / dQ 2 dVL, and the polarized cross section, 
dAa(s T )/dQ 2 dtt. 

The unpolarized Drell-Yan cross section can be factorized according to the QCD factoriza- 
tion theorem |12| , |i~3"|| , 
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±M»£ = g/^ W ^ fl(j0 !^L (6) 

where J2 a ,b run over a ^ possible parton flavors, the 0's are normal parton distributions, and 
do jdQ 2 d£l are perturbatively calculable partonic hard parts. At the lowest order, the partonic 
hard part is given by the qq — > J*(Q) — * ii tree-level diagram, and can be expressed as 

d&aS-+fI(0 > \X 9 ( OL 2 \ „,.„, i is r ,„ ~ „s ( 1 



'qq-*U{Q)X _ 



dQ2dn 4 [^j H^(xp, x'p') L^i, Q-l) l—j 5(Q 2 - (xp + x'p') 2 ) , (7) 

where s = (xp + x'p') 2 ~ xx'S with total cm. energy S = 2p ■ p'. The parton level hadronic 
tensor H^ v is given by 

H^(xp, x'p') = \(xpY{x'p'Y + {xpY{x'p'Y ~ 2{xp) ■ {x'p')g^} (8) 

with the explicit color factor (1/3). The corresponding leptonic tensor is 

L^(£, Q-£)=4 [IpiQ - l) v + (Q- l)pl v -£-Q g^] . (9) 

Both hadronic and leptonic tensors show explicit electromagnetic current conservation: 
Q^H^ U = and Q^L^ V = 0, which will be an issue below when we calculate the twist-3 
contributions to the asymmetry. 

By contracting the hadronic and leptonic tensors, we obtain the lowest order unpolarized 
Drell-Yan cross section in the virtual photon's rest frame, 



dQ 2 dQ 9Q 2 ,, 



4 Jdx q{x) dx' q{x') ^- (- \ [l + cos 2 0] 5{Q 2 - xx'S) 



(10) 



where J2q runs over all quark and antiquark flavors and the factorization scale dependence of 
the parton distributions is suppressed. 

The single transverse-spin asymmetry is a twist-3 effect. In contrast to Eq. (|(J), the Drell- 
Yan single transverse-spin asymmetry has the following generalized factorization formula || 

-^J dx 1 dx a T a/A {x 1 ,X 2l 8 T ) J dxtk/B{x) q , (11) 

where T^ A (xi, x 2 \ St) are twist-3 three-parton correlation functions with a dimension of mass 
and dAa(sT)/dQ 2 dQ are perturbatively calculable partonic parts. The explicit 1/Q factor in 
Eq. (|TTD is a signal of a twist-3 effect and takes care of the mass dimension of the twist-3 
correlation functions. In order to extract the lowest order contribution to the Drell-Yan single 
transverse-spin asymmetry, we need to start with lowest order diagrams with two-fermion and 
one gluon lines from the polarized hadron, as shown in Fig. ^| P, |Tl~| , |l4"f . Because of the virtual 
photon of the Drell-Yan process, contributions from Fig. |2] can be written in general as 

I =(^) W^(Q)L^(£,Q-£) (12) 
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with the leptonic tensor given in Eq. (|9]). After collinear expansion of the parton's momenta, 
and separation of fermion traces |3,11,14], the leading contributions to the hadronic tensor W^ v 
from the diagram in Fig. ^| can be expressed in the following factorized form, 



W^(Q) = J24 J dx l dx 2 T^(x 1 ,x 2 ;s T ) J dx' q{x') H^( Xl ,x 2 ,x';Q) 



(13) 



where J2 q runs over all quark and antiquark flavors. The hadronic matrix element T^(xi,x 2 ; St) 
is given by 



T?(x 1 ,x 2 ;s T ) 



P +d Vl P +d U2 e ix ip +y- e i(x2-xi)p+yz 



27T 2tt 



X {P,S T \lp q (0) 



1 



2p^ 



^{y 2 )%{yt)\p,s T ) . 



(14) 



where A a is the gluon field. The partonic part H£ v in Eq. ( |T3"D is given by the four Feynman 
diagrams in Fig. |3] with the quark lines contracted by -7 • p, and antiquark lines contracted 
by \l ■ p' [p|, pTT] , p^] . As we will demonstrate below, the inclusion of all four diagrams are very 
important for preserving electromagnetic current conservation. The diagrams ||(c) and |3|(d) 
play the same role as the twist-3 contribution of diagram Fig. 1 of Ref. [15] in the calculation 
of Ref. fl8fl. 



We will present a similar derivation in a covariant gauge in the next section. 

In n ■ A(y 2 ) = light-cone gauge, the leading contribution of the gluon field A a comes from 
its transverse components (i.e., a = 1,2) fi~6f , and the corresponding gluon field strength is 
given by 



It is important to note that the T^{x\ } x 2 ] St) in Eq. ( p~4|) is not necessarily a twist-3 matrix 
element, and the hadronic tensor W^ u in Eq. ([U|) could have twist-2, twist-3, and even higher 
twist contributions [14]. In the rest of this section, we show how to extract the twist-3 contri- 
butions to the single transverse-spin asymmetry from the W^ u in Eq. (|T3|) in light-cone gauge. 



F + °( y2 -) = nPd p A°(y^. 
Therefore, in the light-cone gauge, we need to make the following replacement, 



(15) 



(x 2 — a;i)p H 



F +C7 {y 2 



(16) 



in the matrix element T a (xi, x 2 ; st) to convert the gluon field A a (y 2 ) into the gluon field 
strength F +a (y2~) } the same as that in the twist-3 matrix element T^-^xi, x 2 ) in Eq. (f|). 

Because of different behavior of the fermion propagators, we derive the partonic tensor H£ v 
in Eq. ([IB] ) in two steps. We first calculate the contributions from the two diagrams in Fig. |3](a) 
and 0(b). With the Lorentz index a in the transverse direction, we obtain 



H^(a + b)=g 



X 2 — X\ 

x 2 — X\ + it 



5(Q 2 -x 2 x'S)+g» a 



P 



.i" 



6(Q 2 - Xl x'S) 



(17) 



Clearly, the partonic part from the two diagrams violates electromagnetic current conservation. 

After the collinear expansion, the fermion propagators in the other two diagrams in Fig. |3] 
are on mass-shell and contributions from these two diagrams are divergent. As shown in 
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Refs. [p[|14[ , the divergent part of these propagators corresponds to the long-distance contri- 
butions that should be included into the twist- 2 quark distribution; and these propagators also 
have a finite contact contribution, which is one twist higher. The concept of special propagator 
was introduced in Ref. [14]] to identify the high twist piece of each propagator. Graphically, the 
special propagator is represented by a normal propagator line plus a bar as shown in Fig. |3](c) 
and fj(d). The Feynman rule for a special propagator depends on the gauge choice because of 
the nature of high twist contributions. In light-cone gauge, for example, the rule for a special 
quark propagator, merged from incoming quark and gluon lines in Fig. 0(c), is given by 13X 



27 • (x 2 p) 27 • n 

(X2P) 2 + ie 2x2P ■ n 



x 2 - x x 



X2 — x\ + ie 



18) 
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where the ll ie n can be either derived following the example given in the next section, or obtained 
directly from the prescriptions given in the Appendix of Ref 
propagator in Eq. (JlS| 
1(d), 



With the special quark 
we derive the twist-3 contributions from the diagrams in Fig. 01(c) and 



H^(c + d) = g 



x 2 - Xi 
x 2 - Xi + ie 



P_ 

x 2 , 



5(Q 2 -x 2 x'S)+g» 




5(Q 2 - xix'S) 



(19) 



Combining the short-distance twist-3 contributions from all four diagrams in Fig. 0, we obtain 



X2 — X\ 

X2 — X\ + ie 



&[-,--] S(Q 2 -x 2 x'S) 
\x' X 2 J 

+ 9* {^-^X 5{Q 2 - Xl x'S) 

\X' X\ J 



(20) 



The partonic part H£ v with complete twist-3 contributions from all four diagrams in Fig. |3] 
clearly respects the electromagnetic current conservation. 

Having the partonic part H£ v in Eq. (|20| ) and the matrix element T CT in Eq. flU]) with its 
gluon field replaced by the corresponding field strength as in Eq. (0), we derive the twist-3 
contributions to the hadronic tensor W iu (Q) in Eq. (13), 



W^(Q) = Y, e l J dxT ( / q \x,x) J dx'q{x') \ 5{Q 2 - xx' S) | e STann 



x 



afP P'Y . u(P P'Y 



In deriving the above result, we took the imaginary part of the single pole, 

1 



X2 — x\ + ie 



-7riS(x2 — Xi) 



[21) 



(22) 



such that the i from the imaginary part cancels the i in Eq. ([16]) to ensure a real contribution to 
the single transverse-spin asymmetry || . The apparent factor x 2 — x\ in numerator of Eq. ( PD| ) 
cancels the 1/ [x 2 — X\) in Eq. ([IB]) , which is the result of converting the matrix element T a to 
the twist-3 correlation function Tp^\ It is clear that the leading order twist-3 hadronic tensor 
W^ V (Q) in Eq. ([21]) respects electromagnetic current conservation, 
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W^{Q) Q I > = W" V (Q) Q u = 



(23) 



with the virtual photon's momentum = (xp)^ + [x'p 

By contracting the hadronic tensor W flu (Q) with the leptonic tensor L^ v 



in Eq. 



and 



using Eq. 
Fig- @, 



(Jl2p, we derive the leading order twist-3 contributions from the general diagram in 



/(3) = (t4) E e «/ dxT%\x,x) J dx'q(x') 



- jr SQ\sm 29 sin 0)<5(Q 2 - xx'S) 



(24) 



where YJ„ runs over all quark and antiquark flavors. Multiplying the twist-3 contributions in 
Eq. (p4j) by the overall flux, coupling, and phase space factor, (a 2 em /AS) = 1/(25) e 4 (1/327T 2 ), 
as well as a color factor (1/3), we obtain the leading order polarized cross section, 



dAa AB ^u {Q)x (s T ) _ Anal 



dQ 2 dQ 



9Q 2 
1 



x 



Air 



q 
3 
4 



dx Tp (x, x) dx' q(x') 
— j- sin 29 sin <5(Q 2 — xa/S 1 ) 



(25) 



Use T q (x, x) = —T Fi 



(V), 



x,x) 



and divide the polarized cross section in Eq. (|25|) by the unpolarized 



cross section in Eq. (|10[ ), we derive our result for the single transverse-spin asymmetry for Drell- 
Yan lepton angular distribution, 

±E g e 2 q fdxT q (x,x)q(Q 2 /xS) 

Q 



A 



N 



Va 



Tca s 



sin 29 sin <fr 
1 + cos 2 9 



T e 2 



J dxq(x) q(Q 2 /xS) 



(26) 



where Y^ q runs over all quark and antiquark flavors. If we suppress the sum over quark flavors, 

as what was done in Refs. J^J§, the asymmetry A N in Eq. ( p6|) is reduced to the A^f MT ^ in 
Eq. (|). 

With three partons from the polarized hadron, the general Feynman diagram in Fig. [| can 
have both soft-gluon as well as soft-quark poles 0. Because of the infrared behavior of soft 
gluons, the soft-gluon poles from the general diagram in Fig. § can be either a single pole, 
like the l/(x 2 — xi + ie) in Eq. (pDp, or a double pole, such as 1/(22 — x\ + ie) 2 . The single 
pole leads to a dependence on T(x, x) while the double pole results into the derivative term 
xdT(x, x)/dx. For example, consider a diagram in Fig. [|, which gives a higher order correction 
to the Drell-Yan single transverse-spin asymmetry if the photon is virtual. If the photon is 
real, it gives a leading order contribution to the single transverse-spin asymmetry in prompt 
photon production 0. In light-cone gauge, this diagram potentially has a double soft-gluon 
pole: one from the gluon propagator and the other from the gluon attached to the polarized 
hadron. However, in general, one of the two potential poles might be canceled by a vanishing 
numerator. But, if the outgoing gluon in Fig. ^ has nonvanish transverse momentum, both soft 
poles could survive and lead to a double soft-gluon pole 0. 

However, the potential double poles never survive at the lowest order single transverse-spin 
asymmetry in the Drell-Yan angular distribution calculated in this paper. Consider the diagram 
in Fig. 0(a), we have potential poles from the quark propagator, 

. (x\ — £2)7 " V — x 'l " p' 
—(xi — X2)x'S + ie 
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and the i/(x2 —Xi)p + in Eq. (|T6|) from the gluon attached to the polarized hadron in light-cone 
gauge. In terms of the generalized factorization theorem |3,11], momenta of all partons entering 
the partonic short- distance hard parts have only components collinear to their respective hadron 
momenta. For the leading order Drell-Yan process, the incoming quark and antiquark lines in 
Fig. |3] have to be contracted with 57 ■ p and -7 ■ p' to ensure the twist-3 contributions and 
color gauge invariance at this twist. Any other combination will result into power suppressed 
corrections [TJj] . Because the index a in Fig. is transverse in light-cone gauge, the 57 ■ p' 
from the twist- 2 antiquark line eliminates the 7 ■ p' term of the quark propagator, and leaves 
the partonic contribution proportional to (X2 — X\)/{x2 — X\ + ie). The vanishing numerator 
factor (x 2 — X\) cancels the potential 1/(2:2 — ^1) pole from the incoming gluon, and leaves a 
single soft-gluon pole. The derivative term in Eq. @ was obtained because the authors did 
not contract the incoming quark and antiquark lines with the -7 • p and ^7 • p' . For a detailed 
discussion on this issue from the perspective of gauge vector (n) independence, see Ref. |§. 



III. ASYMMETRY DERIVED IN A CO VARIANT GAUGE 

In deriving the single transverse-spin asymmetry in the Drell-Yan angular distribution in 
the last section, we kept explicit electromagnetic current conservation. In order to test the 
color gauge invariance of our result, we present a derivation of the same asymmetry in a color 
covariant gauge in this section. 

Similar to our calculation in light-cone gauge, we start with the general decomposition 
in Eq. fll2|) , and a complete set of partonic Feynman diagrams at the order of 0(ga 2 em ) in 
Fig. [3[ Due to the change of the gauge, one difference is that the gluon field A a (y2~) in the 
matrix element T°(xi, X2', sr) is no longer dominated by its transverse components. Instead, 
the leading contribution of the gluon field is from the component parallel to the direction of its 
momentum, A° oc p a fl6|| . Therefore, we can get the leading contribution from the gluon field 



in a covariant gauge by the following replacement, 

p+ 

In order to convert the A + (y2~) to the corresponding field strength F +P (y2~), we need to give 
the gluon field a small transverse momentum fey, and then convert k T A + {y 2 ) into F +P (y2~) 

With the gluon carrying a small transverse momentum fey, the hadronic tensor W ,J,V {Q) in 
Eq. ( |12"D has a slightly different expression, 

W^iQ) = 114 J dx i dx 2 J d 2 k T T q (x 1} x 2 , k T ; s T ) J dx' q{x') H pu (x u x 2 , k T , x'; Q) , (28) 

where J2q runs over all quark and antiquark flavors. The hadronic matrix element 
T q (xi, x 2 , k T ; st) is defined as 

T q ( Xl , x 2 , k- s T ) = I ^ ^ e^-*^ e** 



x (p,s T \ip q (0,0 
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T 



2p^ 



-^A + (y^,y T ) 



MVi)\p^t). (29) 



S 



The partonic part H^ v {xi, x 2 , kr, x'; Q) is again given by the four Feynman diagrams in Fig. [3], 
except the incoming gluon momentum is now [x 2 — x-y)p + h?, and the gluon polarization a is 
contracted with p° '. 

In order to extract the twist-3 contributions from these diagrams, we need to get one 
power of kx from the partonic part to convert the gluon field A + in the matrix element into 
corresponding gluon field strength ||. Expand the partonic part at kx = 0, 



« H^(x u x 2 , k T = 0, s T ) + -J^H» v (xx, x 2 , k T = 0, x'; s T ) K + Oik, 



(30) 



The first term in the above expansion corresponds to the eikonal line contribution to the twist- 
2 quark distribution. The second term in Eq. (^) is responsible for the twist-3 contribution, 
while the remaining terms correspond to even higher twist contributions. 

Substituting the second term of the above expansion into Eq. (0), and integrating over 
/ d 2 kT, we obtain the twist-3 contributions to the hadronic tensor in a covariant gauge, 



dx x dx 2 T^{ 



Xi, x 2 ) \ dx' q(x') 



where the twist-3 correlation function T^\x%, x 2 ) is given in Eq. (|J). Therefore, calculating 
the twist-3 short- distance contribution is equivalent to extracting the terms linear in k^ from 
the Feynman diagrams in Fig. (|3]) with incoming gluon momentum (x 2 — x i)p + kT- 

As in the previous section, we calculate the short-distance partonic contributions to the 
asymmetry in two steps. We first evaluate the contributions from the diagrams in Fig. |3](a) 
and H(b), and find 



psxnn 



2vr 



d_ 



k T =0 



(31) 



d 



(a + b) = -g 



1 



x 2 — X\ + ie 



#(-,) S(Q 2 -x 2 x'S) 



+g u p ^ys(Q 2 -x 1 x'S) 



(32) 



Again, the contributions from these two diagrams do not respect electromagnetic current con- 
servation, and are exactly the same as what we calculated in the light-cone gauge. Therefore, 
we need to find the other half of the contributions from the diagrams in Fig. |3|(c) and 0(d). 

In a covariant gauge, the gluon polarization index a in Fig. [3] is contracted with p° ' . For 
the two diagrams in Fig. [|(c) and |3](d), the gluon interaction vertices are directly connected to 
the incoming quark lines, which are to be contracted with |7 • p. Naively, we would conclude 
that contributions from these two diagrams vanish because p a [y CT \^ ■ p) oc p 2 = 0. But, due to 
the vanishing denominator of the quark propagators of these two diagrams, we actually have a 
"0/0" situation for the contributions of these two diagrams. 

Same as in the light-cone gauge, these two diagrams have not only twist-3 contributions but 
also long-distance twist-2 contributions. In order to extract the correct twist-3 contributions 
from these two diagrams in covariant gauge, consider a simple and generic diagram in Fig. [5], 
where the incoming quark of momentum X\p is on mass-shell, and the gluon's momentum 
k — > (x 2 — x\)p + kT- The contribution from this generic diagram can be written as 



H = Tr 



. N Z7 • (xip + k) 

Mixip + k) — y~ , (-17,7)7 • (xip) 

[x\p + ky + ie 



(33) 
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where M represents the blob in Fig. ^. Since we are interested in the twist-3 contribution in a 
covariant gauge, we need only the terms linearly proportional to h?, 



H — ► # (3) = Tr [M(x 2 p)7 ■ k T ] 



X\ p 



(x 2 - xi)(2xip 2 ) + ie 



1 



Tr [M(x 2 ph ■ k T ] - 



1 



2 x 2 — X\ + ie 



(34) 



In deriving the first line in Eq. (0), we dropped a k\ in the denominator and the /c^-dependence 
in M, because they correspond to twist-4 or higher contributions. 

Following this approach, we derive the twist-3 contributions from the two diagrams in 
Fig. |(c) and |(d), 



d 



{c + d) 



1 



x 2 — x\ + ie 



P_ 

x 2/ 



5(Q 2 - x 2 x'S) 



8(Q 2 - Xl x'S) 



(35) 



By adding all twist-3 contributions of the four diagrams in Fig. [3], we obtain in a covariant 
gauge 



d 



1 



x 2 — x\ + ie 



9p [-,--) 6(Q 2 -x 2 x'S) 
\x' x 2 J 



x 



(36) 



which clearly respects electromagnetic current conservation. 

By substituting Eq. (^) into Eq. (|3~1"D, and contracting with the leptonic tensor L^ u , we 
derive the same polarized Drell-Yan cross section in Eq. (|25|) and the same asymmetry in 
Eq. ( Pop in a color covariant gauge. Therefore, we can conclude that our result in Eq. (^) 
for the single transverse-spin asymmetry in the Drell-Yan lepton angular distribution not only 
respects electromagnetic gauge invariance, but also is the same in both color light-cone and 
color covariant gauges. 



IV. NUMERICAL ESTIMATE AND CONCLUSIONS 

In order to estimate the numerical size of the single transverse-spin asymmetry in the Drell- 
Yan lepton angular distribution in Eq. (|26|), we need to know the size and sign of the twist-3 
quark-gluon correlation functions T^\x,x). 

Fortunately, the same twist-3 quark-gluon correlation functions T^\x, x) appear in a num- 
ber of other single transverse-spin asymmetries |||6| . In principle, self-consistency of the twist-3 
correlation functions in different physical observables is a test of QCD and its factorization the- 
orems beyond the leading twist physics. 

In Ref. 0, single transverse-spin asymmetries in hadronic pion production were calculated 
in perturbative QCD, in terms of the same generalized factorization theorem that was used 
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in this paper. The same twist-3 quark-gluon correlation functions Tpf\x,x) dominate the 
leading contributions to the asymmetries. Although the existing data on pion asymmetries 
have relatively low pion transverse momenta, the general features of the data are naturally 
explained by the theoretical formula's ||. By comparing the theoretical calculations with 
the existing data on the difference between 7r + and tc~ asymmetries, it seems that the twist- 
3 correlation function Tpj (x, x) for valence up quark has an opposite sign to T^P (x, x) for 

valence down quark. In addition, a positive Tp^\x,x) is favored. 

If we adopt the model introduced in Ref. || for the twist-3 quark-gluon correlation functions, 

T%\x,x)*K q \q(x) (37) 

with k u = +1 = —Kd, k s = 0, and A ~ 100 MeV, we estimate the single transverse-spin 
asymmetry in the Drell-Yan lepton angular distribution as follows. Taking the angle ~ 90 
degrees, we have the angular dependence, 



sin 28 sin . 



1 + cos 2 9 



< 0.7. (38) 



Because the valence up quark is about twice the valence down quark, and the fractional charge 
square e 2 = 4/9 is a factor of 4 larger than e 2 d = 1/9, we approximate 

J2 q e 2 q JdxT q (x,x)q{Q 2 /xS) x fdxu{x) u{Q 2 /xS) 



A 



E q e 2 . / dx q{x) q{Q 2 /xS) j dx \u{x) u(Q 2 /xS) + u{x) u{Q 2 /xS)} 

where the minus sign is due to the definition T q (x, :). Because of the assumption, 

k a = 0, we have only one term in the numerator in Eq. (^) in comparison with the two 
terms in the denominator. The 2 in the second line of Eq. fl3"P| ) is due to the fact that the two 
terms in the denominator are equal for proton-proton collision. Taking the maximum value 
from the angular dependence in Eq. (|3"8"D, we estimate that the magnitude of the asymmetry in 
Eq. (p6|) can be as large as (for y/4na s ~ 2) 

|^|~0.7^. (40) 

For A ~ 100 MeV, we have A N of about 3.5% at Q = 2 GeV, and 1.75% at Q = 4 GeV, where 
we chose the Q values to be just below and above the J/ip resonance. Considering that RHIC's 



sensitivity to double transverse-spin asymmetries in Drell-Yan is at the percent level [18], these 
values for the single transverse-spin asymmetry might be detectable. 

The numerical values for the single transverse-spin asymmetry in Drell-Yan lepton angular 
distribution are small, which is consistent with the fact that the single transverse-spin asym- 
metries are a twist-3 effect. For any inclusive observables, like the Drell-Yan cross section 
discussed here, the asymmetries should be of the order of Aqcd/Q with Q being the energy 
exchange of the collisions. However, the single transverse-spin asymmetries can be large for 
observables at certain phase space where the twist-2 contributions are steeply falling, and the 
twist-3 contributions get a kinematic enhancement due to the derivative term 0. 
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We would like to contrast our small asymmetry estimates with those of Ref. Q, which find 
large asymmetry values based on their model and its parameters fitted to the pion production 
asymmetries. A future measurement of the size of the asymmetry can therefore clearly distin- 
guish between the model predictions of [|J (asymmetries between 20% and 40% for the whole 
xf range) and the estimates presented here. 

In conclusion, we have calculated the single transverse-spin asymmetry for the Drell-Yan 
lepton-pair's angular distribution in perturbative QCD. We found that our result at leading 
order in the strong coupling constant, given in Eq. fl26|) , is consistent with what was derived 
in Ref. ||, but different from what was found in Ref. (J7|. We derived our result in both color 
light-cone and covariant gauges while keeping explicit electromagnetic current conservation. 
With a model of the twist-3 correlation function, our calculated asymmetry predicts both the 
sign and magnitude of the asymmetry. 
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FIG. 1. Coordinate choice for the virtual photon's rest frame. 




FIG. 2. General diagram contributing to the Drell-Yan single transverse-spin asymmetry. 
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FIG. 3. Lowest order Feynman diagrams contributing to the Drell-Yan single transverse-spin 
asymmetry. 
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FIG. 4. Sample Feynman diagram generating a double soft-gluon pole. 




FIG. 5. A generic diagram for extracting the special quark propagator. 
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